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I present a solution to the full Einstein-fluid equations representing a self- 

q'. gravitating Bjorken flow. The motion and the geometry become inhomogeneous 

in the plane transversal to the flow and the energy density profile acquires, due 
to gravity, corrections in terms of proper time as compared to the original test 
hydrodynamics. The transverse distribution of energy density, for example, 
becomes e(r, r)/e(r, 0) = cosh - (3ar). 
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Bjorken flow [lj represents the most fascinating application of relativistic hydro- 
dynamics to an extremely complex physical system describing an average motion of 
partons resulting in a collision of heavy ions. The application of hydrodynamics to 
similar problems was pioneered by Landau [2] to describe the high-energy multipar- 
ticle collisions. Both in Bjorken and Landau descriptions it is assumed that after the 
collision of heavy ions the mean free path of the constituencies is short enough, so 
that the hydrodynamical description is meaningful. The difference between the two 
pictures is in the symmetry assumptions. In Bjorken hydrodynamics, one assumes the 
so-called boost invariance, so that the energy density only depends on proper time r, 
while in Landau picture, no such symmetry restriction is made and the density may be 
a function of all spatial coordinates. Needless to say that Bjorken flow is a particular 
and simpler version of Landau hydrodynamics, nevertheless, it is surprising it works 
so good [3J. A different and renewed motivation in these studies comes from their 
relation to the AdS-CFT correspondence conjecture [I], [5], [6], because they serve as 
an input to understand the highly nontrivial behavior of Quantum Chromodynamics 
in a strong coupling regime. 

The symmetries one imposes on the Bjorken flow are: the boost symmetry along 
the beam, translational and the rotational invariance in the transverse plane. These 
symmetries allow one to parametrise the "future wedge" of the Minkowski spacetime 
in the following way 

ds 2 = -dr 2 + r 2 dr] 2 + dp 2 + p 2 d<\> 2 (1) 

Here r is the proper time and 7] is historically called the rapidity and the rest are 
usual cylindrical coordinates. 

For further purposes I will write the metric in the following form 

ds 2 = -dr 2 + dp 2 + rp (r/p dr] 2 + p/r dcj) 2 ) (2) 

The spacelike Killing fields Jr and ■£- are the rotational and the boost Killing 
vectors respectively, and the form of the line element ([2]) gives one an idea as to how 
to proceed in order to generalise this set up to a General Relativistic flow. 

Now, if one introduces an ideal fluid with the linear equation of state p = 1/3 e 
with the energy density and the pressure such that these only depend on the proper 
time r, in other words the fluid velocity has no tilt, it is quite easy to solve the 
hydrodynamical equations and to obtain that the fluid density scales as oc r _4//3 . 
This is quite close to the experimental picture, grosso modo, within the range where 
the hydrodynamics makes sense. 

On the other hand, the above picture is quite idealised, of course, some amount of 
viscosity should be added and some symmetries relaxed [7], nevertheless, it is quite 
surprising that the picture works so well. 

While the Bjorken hydrodynamics deals exclusively with the flow on a given fixed 
flat background geometry (test hydrodynamics), and given that even slight relaxations 



of symmetry would leed to quite complicated nonlinear hydrodynamical equations, 
even the so-called Khalatnikov solution [8] of a 1+1 dimensional flow is quite a mess 
[9], it is yet another pleasant surprise that one may integrate an exact general rela- 
tivistic solution which describes self gravitating Bjorken flow. The main purpose of 
this Letter is to present such a solution and to compare it to the original Bjorken's 
test hydrodynamics. As a by-product, I will also obtain solutions to the test hydro- 
dynamics without tilt on a class of cylindrical geometries. 

I will stick to the original Bjorken picture as close as possible and will assume that 
the fluid velocity has no tilt (see, however, a comment after the equation (jfJJ).) Never- 
theless, since one must solve selfconsitently the coupled Einstein-fluid equations, one 
can not expect that the geometry would share all the above mentioned symmetries. 
This is the essence of the Einstein theory, the matter influences the geometry which 
in turn changes its motion. We may keep the boost and the rotational Killing vectors 
intact, however, there is no reason why the line element should not depend both on 
t and p coordinates. In fact, the form of the line element ([2]) indeed suggests the 
dependence on p. I will therefore assume the following geometry 

ds 2 = / ( -dt 2 + dr 2 ) + g (q dif + q~ l d(j) 2 ) (3) 

Here /, g and q are functions of both t and r the "conformal" coordinates which 
are labled differently to distinguish them from the proper time and the proper distance 
coordinates. Obviously for the original Bjorken flow these functions are / = 1, g = tr 
and q = t/r. The equation ([3]) seems a natural generalization of Bjorken geometry, 
while the attempts to use the homogeneous Kasner, or flat FRW line element, as 
some authors do, fails to address the inherited symmetries of the problem. 

I now specify the matter. The perfect fluid is assumed to have a linear equation 
of state, and because the fluid flow is irrotational, which is a must in this geometry, 
one may introduce the following velocity potential a |10| . [TT] : 
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As will bee later seen, the velocity potential is an extremely useful tool to solve 
the hydrodynamics. 

Having done so, one may further define the kinetic scalar ("enthalpy") X = 
— 1/2 a a a a . The pressure and the energy density can then be expressed as follows 



p = p(X), e = 2Xp'-p (5) 

Here prime, as usual, stands for the derivative of the function with respect to its 
argument. 

If the equation of state is linear p = we, one may further write 111] . [T2] 



p(X) = X {w+1)/2w (6) 



We now assume that we have chosen our coordinates comoving with the fluid flow 
as in the original setting so that the velocity has only a zero component uq. This 
translates, in terms of the velocity potential, that a is a function of t alone. In fact, if 
one would have even allowed a tilt, so that the velocity would "catch" a component in 
the transversal direction (§^7^ 0), as for example in (|7j), there would still be a way 
to introduce a new coordinate system comoving with the fluid, and maintain the form 
of the metric [T3]. This would not necessarily be true for the test hydrodynamics, 
and it may also spoil the separability of the metric functions which I will assume in 
the future. 

The full Einstein equations for the line element ([3]) with the fluid specified above 
are found in [11]. It is instructive, however, to display the dynamical equation for the 
velocity potential V^ (p'cr M ) = [12]. This reads: 

1 1 / 1 \ 

& = 0, (7) 

where the velocity of sound v s is given by 
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and can be easily obtained from the relations (5) and the expression v 2 = -^-. 

Note that the transversal degree of freedom of the metric q plays no explicit role 
in the dynamics of the fluid, but does influence the flow via the full Einstein equa- 
tions contributing to the longitudinal expansion /. Both, though, the longitudinal 
expansion / and the function g, which is proportional to the area of the isometry 
group orbits, do appear in the equation. These functions (p, g and /) are determined 
by the Einstein Equations. 

Using the linear equation of state (jH]), so that v 2 = w, and assuming that all the 
functions of the metric are separable, f — fr fa, 9 = 9t 9r and so on, one may 
easily integrate the dynamical equation ((Tj) to get 

- = b^~, (9) 

9t 

where lower index T indicates the time dependent part of the respective function 
and b is an arbitrary integration constant. The solution §§§ represents therefore 
Bjorken flow on a generalised geometry given by the line element ([3]). Hence, given 
the background geometry Q with separable functions, the velocity potential is given 
by the solution (Q, and therefore all the kinematical fluid variables, may be easily 
evaluated. For the Bjorken geometry, and w = 1/3 we get a oc t 2 ' 3 which leaves us 
with e oc £~ 4//3 and t is a proper time coordinate in this case. For the spatially flat 
FRW geometry sourced by the same fluid (radiation) / oc g oc t 2 , we get a oc t, 



X oc t~ 2 and e oc t -4 which becomes e oc r -2 in terms of proper time. One may play 
round with the equation (J7J) further, but I will leave it for elsewhere. 

The equation (|7|), though important, is only a part of the full Einstein Equations 
and describes a test hydrodynamics on a given geometry. If one solves the rest of 
the equations with the following initial and boundary conditions compatible with the 
Bjorken geometry, g — > tr as t and r — > 0, and t — > as r — > one can integrate 
the following solution to the Einstein Equations: 

/ = sinh 4 (at) cosh 2 (3ar) 

g = sinh (at) sinh (3ar) cosh -2 ' 3 (3ar) 

g = sinh (at) sinh (3ar) (10) 

The velocity potential for this solution is given by 

a= (15 a 2 ) 1/4 cosh (at) (11) 

while the energy density becomes 

e = 15a 2 sinh" 4 (at) cosh" 4 (3ar) (12) 

Here a is a free parameter specifying the density. In fact this solution was derived 
some years ago by the author together with J. Senovilla [2] in the context of inho- 
mogeneous cosmology, but it could not occur to us then that the solution describes 
the inhomogeneous self gravitating Bjorken flow. Indeed, since the coordinate t is not 
the proper time, we find that the coordinate time t oc r 1 / 3 and therefore the energy 
density scales as 

e oc r" 4 / 3 , (13) 

to the lowest order in proper time. However, some corrections are due. Because 
we have chosen our coordinates in a way that r remains constant along the fluid lines, 
the proper time is given by the following expression 

r= fffdt=±- (at) 3 + -L (at) 5 + 0((at) 7 ) (14) 

jo oa 15a 

The energy density e then evolves as (just the first two terms) 

e oc a (a) r^ 13 - (3 (a) r -2 / 3 , (15) 

here a and j3 are both positive functions of the parameter a, I have assumed a constant 
r, and have used the lowest order of the proper time expansion. Of course, having at 
hand the exact solution, there is no need in the series expansions, nevertheless, these 
are instructive in order to further elucidate the physics. As one can easily see from the 
equation f[T5l) . the second term becomes dominating at late times, where the energy 



density appears to become negative. This is an artifact of the series expansions, for 
the exact energy density never becomes negative. On the other hand, while the first 
term scales as the energy density of the ideal fluid, the second term acts as if it were 
viscosity by enhancing the fall off of the evolving density. Of course, the fluid remains 
inviscid for all the time, and this effect is purely due to self- gravity. Assuming the 
Stefan-Boltzmann's law (e ~ T 4 ) one can easily find the temperature distribution 
and then define the temperature contrast ST as 

T-T R 
ST = —A (16) 

J-B 

where T is the temperature found from the exact solution, while Tr represents 
the temperature of the test flow. This temperature contrast evolves as ST ~ — t 1 / 3 . 
Another interesting physical quantity is the distribution of energy density in the 
transverse coordinate 

e(r,r)/e(r,0) = cosh" 4 (3ar), (17) 

which is given by a neat simple expression in terms of "conformal" distance r. It 
is assumed that this quantity is proportional to the distribution of nucleous in the 
fireball, the actual numbers, of course, would depend on the colliding constituencies. 

To close, I have presented an exact solution to the Einstein-fluid equations which 
describes selfgravitating Bjorken flow. The gravity changes the energy density distri- 
bution and its evolution in proper time. From the exact expressions one may easily 
find all relevant kinematical and thermodynamical quantities. As a by-product, I 
have obtained some test hydrodynamical solutions on the expanding cylindrical ge- 
ometry. It would be interesting in the future to study more the test hydrodynamics 
with nonlinear equations of state, as well as to consider the solution ffTUj) as an in- 
put to study the ADS/CFT correspondence in situations where the geometry is both 
inhomogeneous and evolving in time. 
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